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We investigate chiral symmetry restoration in finite spatial volume and at finite 
temperature by calculating the dependence of the chiral phase transition temperature 
T c on the size of the spatial volume and the current-quark mass for the quark-meson 
model, using the proper-time Renormalization Group approach. We find that the 
critical temperature is weakly dependent on the size of the spatial volume for large 
current-quark masses, but depends strongly on it for small current-quark masses. In 
addition, for small volumes we observe a dependence on the choice of quark boundary 
conditions. 

I. INTRODUCTION 


Phase transitions in Quantum Chromodynamics (QCD) are currently very actively re¬ 
searched. Most of the attention is focused on the phase transition at finite baryon density 


and temperature 




where the existence of a critical point in the phase 
. Even at vanishing density, the order of the 


diagram is not yet conclusively settled [4,J2, 
phase transition is still under discussion [llj, 

While QCD is perturbative at large momentum scales, the low-energy limit of the the¬ 
ory is dominated by non-perturbative phenomena. This makes non-perturbative methods 
indispensable, in particular for investigations of t he p hase transition. Effective low-energy 
theories such as chiral perturbation theory 


3ns of t he p i 

B,HQc 


describe the low-energy limit of QCD 
well, but cannot address the restoration of chiral symmetry and the deconfinement tran¬ 
sition: These questions require a connection to the high-momentum degrees of freedom. 
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Lattice gauge theory, on the other hand, yields non-perturbative results, allows an investi¬ 
gation of the phase transition, and can in principle provide the necessary effective couplings 

But even in the light of recent 


, most lattice simulations still require 


for the description by a low-energy effective theory 
advances with light fermions on the lattice |2y, |2l|, 
extrapolations to small, realistic pion masses. In any case extrapolations to the continuum 
limit and to infinite volume are necessary. In particular, finite volume effects are more severe 


when the pion mass approaches the chiral limit 


22 


23 


m 


M 


Therefore, the influence 


of the finite volume for small pion masses should be studied with other methods in parallel. 

In addition, lattice gauge theory provides little guidance to understand the emergence of 
low-energy dynamics. The interplay between lattice gauge theory and other non-perturbative 


met hoc 
Group 


s such as Dyson-Schwinger equations [23, I2& 12-.1 


ill and the Renormalization 


32 


33 


3)4 i35l I -‘I fil l should prove fruitful to further our understanding. There is also a 


need for model systems that describe particular aspects of the dynamical generation of the 
low-momentum physics from the high-momentum theory. One example is the description 


0 


of chiral symmetry breaking via the Nambu-Jona-Lasinio model laa and its modifications. 
The quark-meson model that we use in the present work belongs to this class of models U). 

Of course, such a model approach cannot answer questions outside the applicability of 
the model. For example, the order of the phase transition is in our case already fixed 
by the 0(4)-symmetry of the model, while the order of the transition in QCD has not 


yet been unambiguously determined 


We must also limit our investigation to the 


chiral phase transition. In QCD, there is no requirement that the chiral and deconfinement 
phase transitions occur at the same point [^j], although so far there is no indication for 
two transitions. On the other hand, the model has been successfully used to investigate 


the quark mass dependence of the chiral transition m and the critical behavior at finite 


density py with Renormalization Group methods. It has also recently been combined with 
Polyakov loop results from the lattice to describe thermodynamical observables from lattice 
QCD [40]. 

The Renormalization Group (RG) is an important tool for the investigation of non- 


perturbative physics |41|, 


42 


m 


44j. In particular, RG methods are well suited to describe 


physics across different momentum scales, and generation of the low-energy effective theories 
from the dynamics. While we do not directly address this issue here, the study of critical 
behavior is of course also well within the scope of an RG approach, and in the present context 
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it has 
139 


reen applied to determine critical exponents, for example for the quark-meson model 


45]. 


In this paper, we consider the chiral phase transition in the framework of the quark- 
meson model. We will apply Renormalization Group methods to calculate the transition 
temperature, its dependence on the quark mass, and its dependence on the size of a finite 
volume. In addition, we investigate the effects of different boundary conditions for the quark 
fields. In calculations based on an effective field theory like chiral perturbation theory, chiral 
symmetry breaking is assumed from the beginning and the values of the effective low-energy 
constants are fixed. In contrast, in our model chiral symmetry is broken dynamically and 
effects of the finite volume on quark condensation are taken into account. We believe that 
such effects could still be important in simulations at the current lattice sizes of order 
L ~ 2 fin, in particular for realistic quark masses. 

In finite volume, strictly speaking no phase transition is possible, since non-analyticities 
cannot appear in the thermodynamic potential (see e.g. |4f|). In general, the investigation 
of phase transitions and critical behavior from results obtained in a finite volume is difficult 
and requires an extrapolation to the large-volume limit. In addition, if a symmetry is 
restored across the transition, this usually requires the introduction of an external held which 
explicitly breaks the symmetry. Even if there is no true order parameter that vanishes strictly 
in one of the phases, rapid changes over a small temperature range are an indication of a 
(crossover) transition. Often, peaks in susceptibilities or other higher-order derivatives of the 
thermodynamic potential are used as criterion to define a pseudo-transition. Here we propose 
to use the mass of the scalar mode, which corresponds to the inverse correlation length for 
fluctuations in the quark condensate, to identify the transition point: A distinct minimum 
of the mass appears at almost the same temperature at which the chiral quark condensate 
drops rapidly. We stress that the implementation of an explicitly chiral-symmetry-breaking 
term is essential, since the finite-volume system will otherwise always be in a regime with 
restored chiral symmetry, once all quantum fluctuations are taken into account. 

Lattice simulations are affected by similar problems. Thus, finite volume effects are 
actually of profound importance in lattice determinations of the order of the phase transition. 


For the determination of the universality class and the critica 
a scaling analysis of thermodynamic observables is necessary 


exponents of the transition, 
L 12L 47, 48L 49]. It remains 


difficult to assess whether current lattice sizes are sufficiently large to observe the expected 
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scaling 


3QQ 


On the other hand, the volume dependence of the (pseudo-) critical 


behavior can be turned into a tool for the analysis: Finite-size scaling of the results has been 


used to test the compatibility of critical exponents with lattice data Ihl Q. We expect 
that future progress in RG analysis will provide much useful insight into these questions. 

The paper is divided into five sections. In the next section, we will give a short review of 
the quark-meson model and the derivation of RG flow equations in the proper-time formu¬ 
lation of the RG. In section mu we will concentrate on the chiral phase transition in infinite 
volume, and its dependence on the parameter of explicit chiral symmetry breaking. In sec¬ 
tion m we will then look at the finite volume effects and the influence of the additional 
momentum scale introduced by the finite volume. We close with a summary and conclusions 
in section El 


II. RG-FLOW EQUATIONS FOR THE QUARK-MESON MODEL 


To determine the chiral phase transition temperature for finite volumes and finite current 
quark masses, we use the chiral quark-meson model 1 . This model is an 0(4)-invariant 
linear cr-model with Nj = 4 mesonic degrees of freedom (a, n) coupled to Nf = 2 flavors 
of constituent quarks in an SU(2)l x SU(2)r invariant way. It is an effective low-energy 
model for dynamical spontaneous chiral symmetry breaking at intermediate scales of k < 
A uv « 1.5 GeV, but it does not contain gluonic degrees of freedom and is not confining. 
The ultraviolet (UV) scale A uy & 1.5 GeV is determined by the validity of a hadronic 
representation of QCD. At the scale A uv, the quark-meson model is defined by the bare 
effective action 

r At „[$]= f d 4 xjg(^ + gm c )q + gq(a + ir • + i(<9 M </>) 2 + U Al/y (0)j (1) 

with a current quark mass term gm c which explicitly breaks the chiral symmetry, and with 
(j) T = (a, 7T 1 , 7r 2 , 7T 2 ). The mesonic potential is characterized by two couplings: 

Uauv ( 0 ) = \ m uv4> 2 + \^uv{0 2 ) 2 • ( 2 ) 


1 In the present approach to the phase transition at vanishing baryon density, we do not include vector 
mesons. The role of the vector meson in medium is not yet completely understood 0EH0E30 
An analysis on the lattice suggests that at high temperature the vector coupling is small compared to the 
scalar coupling 0j. 
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In a Gaussian approximation, we can perform the functional integration of the bosonic and 
fermionic fields and obtain the one-loop effective action for the scalar fields </>, 


I’M = r A „ v M - Tr log (r<?M) + ^Tr log (rg'M) (3) 

/o\ /o\ 

where r { B [(f)\ and T}/[</>] are the inverse two-point, functions. We neglect a possible space 
dependence of the expectation value and take the wave-function renormalization and the 
Yukawa-coupling to be constant. Since the traces in Eq. © are infrared (IR) divergent, we 
use the Schwinger proper-time representation of the logarithms and introduce an infrared 
cutoff function 2 f a (rk 2 ), where the variable r denotes Schwinger’s proper time and k is a 
cutoff scale. The derivative of the cutoff function with respect to the scale k is given by 


fc S / “ (Tfc2) = -rPTI) (Ti 


2\a+l ^—rk 2 


( 4 ) 


The inverse two-point functions in Eq. © depend on the second derivatives of the effective 
potential U. By replacing the bare masses and couplings in the inverse two-point, functions 
with the scale-dependent quantities, we obtain the so-called renormalization group improved 
flow equation for the effective potential 14, in infinite volume for zero temperature and finite 
current quark mass jz^: 


k^-U k (a,n 2 ,T 


fc 2(a+1) / 4 N c N f 

^ OC) ~ 16a(a-l)vr 2 j ~~ (k 2 + M 2 {a, W 2 ))— 1 

1 Nj-1 \ 

( k 2 + M 2 (a, 7f 2 )) a_1 ( k 2 + M 2 (cr, 7r 2 )) a_1 J 


( 5 ) 


Integrating the flow equation from the UV scale to k —> 0, we obtain an effective potential 
in which quantum corrections from all scales have been systematically included. 

Since we allow for explicit symmetry breaking, the O (4)-syrrirrietry of the effective poten¬ 
tial is lost. However, it remains 0(3)-symmetric in the pion-subspace, so that the pion-fields 
can only appear in the combination n 2 on the right-hand side. In order to be able to perform 
the Schwinger proper-time integration in infinite volume, we have to choose a > 2 45j. 

The meson masses M a and M n in Eq. © are the eigenvalues of the second-derivative 


matrix of the mesonic potential 


ial, cf. 0 


for an explicit representation, and the constituent 


2 Although the RG flow equations themselves depend on the particular form of the cutoff function, physical 
quantities calculated from the RG flow should not depend on the choice of the cutoff function in the limit 
k -► 0 . 
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quark mass M q is given by 


M q = g 2 [(cr + m c ) 2 + tt 2 } . (6) 

We generalize the renormalization group flow equations to a finite four-dimensional Eu¬ 
clidean volume L 3 x T by replacing the integrals over the momenta in the evaluation of the 
trace in Eq. m by a sum 


fd Pi ... ^ ± .... (7) 

J rii ——oo 

The boundary conditions in the Euclidean time direction are fixed by the statistics of the 
fields. The thermal Matsubara frequencies take the values 


u> no = 2nn 0 T and u no = (2 n 0 + 1 )vrT, 


( 8 ) 


for bosons and for fermions, respectively, where the temperature is denoted by T. However, 
we are free in the choice of boundary conditions for the bosons and fermions in the space 
directions. In the following we use the short-hand notation 

2 47T 2 2 2 47T 2 / 1\ 2 

pp = jY2^ n i and ^ = h + 2 ) (9) 

i=1 i =1 k / 

for the three-momenta in the case of periodic (p) and anti-periodic (ap) boundary conditions. 
We will consider both choices for the quark fields, but employ only periodic boundary 
conditions for mesonic fields. Then the flow equation for finite temperature and finite volume 
reads 


to 

no n 

Nj=4 

_l_ ^ e -H fc2 +<, +Pp+ M iW 2 )) 
i =1 

The sums in Eq. m run from — 00 to + 00 , where the vector n denotes (tt-i, ri 2 , n^). Since we 
have to solve the flow equation numerically, we rewrite it in terms of Jacobi-Elliptic-Theta 
functions: 



fc^-f/ fc (a,7f 2 ,T,L) = 


fc2(a+l) 


T(a + 1) 


d,TT a -4 N c N f e 


~r{k 2 +ul 


-,+Plp,p + M q (o-,7T 2 )) 


k^- k U t (a,if,T,L) = (- iNcN ,Q%] p [a,(k 2 +M^))L\TL) 


(4 vr )a+ 1 r( a + 1) 

Nj=A 

E Oi B) (a,(k 2 +M?(a,n 2 ))L 2 ,TL 


2—1 


( 11 ) 
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We have introduced the auxiliary (dimensionless) functions 


/ oo / v 3 

dss a e~^d ap (st 2 )(^d ap (s)j , 

r°° / \ 3 

0j P(a,u,t) = j ds s a e~^i9 ap (st 2 ) (tf p (s)J , 

r°° / \ 3 

@^(a,u j,t) = / ds s a e~^'d p (st 2 )^t9 p (s)j , 


( 12 ) 

(13) 

(14) 


where d p and d ap are Jacobi-Elliptic-Theta functions defined as 


$ ap (x ) = X e 

n=—oc 

oo 

^0) = X e 


-X7r(n+i) 2 _ ^,-i 


9=1 


= X 2 


E l 7TP 

(—l) 9 x _2 e * 

OO 

^+ 2 E' 


1_ _ 7rg 

x 2 e * 


(15) 

(16) 


3=1 


The representation in terms of these functions accelerates the numerical calculations by a 
factor of about a hundred, compared to the representation used in Refs. J2 H] The first 
representation in Eq. m and m is the standard Matsubara summation of the momenta. 
The second representation on the right hand side in Eq. m and m is obtained by applying 
Poisson’s formula to the first representation. One can use this representation to separate 
the zero-temperature and infinite-volume contribution of the flow equation. Indeed, using 
the approximation i? ap (s) = i? p (s) « s~^ in Eqs. ED, (USD, and (ED yields 


9iP(a,v,t) = 0 { p F) (a,u,t) = e { p B \a,u,t) « 

Z UJ 


(17) 


Inserting this in Eq. ED, we obtain the flow equation © for infinite volume and zero 
temperature. From now on, we use a = 2 for both the infinite-volume and f ini te-volu m e 
calculations. 

The flow equations © and ED are partial differential equations which can be solved by 
using a projection of these flow equations on the following ansatz for the mesonic potential 


26]: 


(18) 


N a i+j<No 

U k (a,TT 2 ) = Y] X a ij(k)(a - a 0 (k)Y(a 2 + TT 2 - a^(k)y 

i= 0 j= 0 

Such a projection results in an infinite set of coupled first-order differential equations for 
the coefficients a,ij(k). In order to solve this set of equations, we limit the sum in Eq. (HK1) 


by choosing N a — 2 


25 


26]. The boundary conditions for the differential equations for 


the coefficients a lJ (k) are determined at the ultraviolet scale k = A jjy. F° r a given current 








quark mass gm c , we determine the initial conditions dij(Auv) for infinite four-dimensional 
Euclidean volume in such a way that we obtain values of m n and f n = cr 0 which are consistent 
with chiral perturbation theory 3]. Consequently, our calculation cannot predict the values 
of m n and f n = cr 0 in infinite volume, but allows to study the behavior of the masses and 
the pion decay constant in a finite four-dimensional Euclidean volume. The parameter sets 
which we have used for the calculations in Sec. HU and ED are listed in Appendix El where 


the relations between the meson masses and the coefficients are summarized as well. 


III. CHIRAL PHASE TRANSITION TEMPERATURE IN INFINITE VOLUME 


In this section, we discuss the dependence of the chiral phase transition temperature in 
infinite volume on the zero-temperature pion mass m i°' 1 = m n {T = 0). In order to define a 
chiral phase transition temperature in the presence of explicit symmetry breaking, we use the 
dependence of the cr-mass on the temperature. We define the phase transition temperature 
T c through the minimum of the cr-mass, 


dm a {T) 


dT 


= 0 and 


d 2 m a {T) 


T=T n 


dT 2 


> 0 . 


(19) 


T—T r 


Alternatively, one can define the phase transition temperature as turning point of the pion 
decay constant as a function of temperature 3 . We have checked that the values for the chiral 
phase transition temperature T c obtained from these two different definitions agree within 
a few percent. For example, in Fig. 0] we compare the normalized cr-mass R[m a \{T) = 
and the normalized pion decay constant R[f n \(T) = y^fy- for = 100 MeV. One observes 
that both definitions for the critical temperature yield for practical purposes the same result. 

In Tab. □ and Fig. [3 we show the chiral phase transition temperature T c obtained in 
this way as a function of the pion mass We find that the transition temperature T c 

depends on the pion mass in the following way, 


Tc(mf) 


Oq + aimf 3 + a 2 {mf > ) 2 + 


( 20 ) 


3 There is no unique definition for the crossover or the pseudo-critical temperature. For example, the critical 
temperature T c can also be defined as the temperature at which reaches half its zero-temperature value 
[ 5 ^. The results for T c obtained with such a definition will in essence agree with our results, as suggested 
by Fig. □ 

4 We do not show lattice results for comparison in this figure since there is no data available for 


mi 0) < 300 MeV in Ref. 


. Q. 
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T [MeV] 


FIG. 1: Normalized sigma mass R[ma](T) = (o) an d normalized pion decay constant R[f n ](T) = 

^ as a function of the temperature, in infinite volume for a pion mass m ® = 100 MeV. 


where the parameters can be determined from a fit to our numerical results as 


a 0 = 149.58 MeV, a x = 0.24258, a 2 = 0.00029 MeV 


-l 


( 21 ) 


a o is then the value for the chiral phase transition temperature in the chiral limit as obtained 
from the fit. A similar relation was also found in lattice simulations jul [vij] with two or 


three quark flavors. The corresponding relation is 
T c (Nf, m PS ) _ T c (Nf, m P s = 0) 


+ h(N f ) 


m 


2 /PS 
PS 


0(m 


PS) i 


( 22 ) 


'a v^ vcr 

where m P s denotes the mass of the pseudoscalar meson, and the string tension a is used to 
set the scale in the lattice calculation. /3 and 5 are the critical exponents of the 0(4)-model 
in three dimensions. The coefficient h(Nf) depends slightly on the number of quark flavors 

Q. 

The analysis in Eq. (E2D assumes that the transition falls into the 0(4) universality class, 
where the ratio of the critical exponents obeys 1/(38 — 0.55. Then, the first order correction 
term is approximately linear, in agreement with our result. On the lattice, however, the 
coefficient of the approximately linear term is about one order of magnitude smaller than the 
result of our calculation. For Nf = 3, a lattice QCD calculation 13] gives h(N f = 3) ^ 0.039. 
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m n [MeV] 

0 

30 

50 

75 

100 

125 

150 

175 

200 

225 

250 

275 

300 

T c [MeV] 

147.6 

157.5 

163.9 

170.5 

178.1 

184.2 

191.8 

200.2 

208.3 

218.1 

228.0 

238.5 

249.3 

R c {m^) 

0 

0.067 

0.104 

0.155 

0.207 

0.250 

0.300 

0.356 

0.411 

0.478 

0.545 

0.616 

0.689 


TABLE I: Dependence of T c and R c (m^) = Tc ^ m ^ on m^\ 


While the exact value for h(Nf = 2) is not given in |l3[, the authors point out that it is of 
the same order of magnitude as the value for Nf = 3. 

As we will see in the next section, it is not possible to explain the smaller value of the 
(approximately) linear term found on the lattice as a Enite volume effect: Since a finite 
volume effect is more severe for smaller pion masses and since it leads to a significantly 


reduced transition temperature in our model, we expect that the slope of T c (mi°' ) ) should 
actually increase in a finite volume, compared to the infinite-volume result. We think that 
the discrepancy may be a consequence of neglecting the gauge sector in the quark-meson 
model. In the chiral limit, the chiral phase transition temperature on the lattice is about 


30 MeV larger [3] than the value obtained in the quark-meson model. 

Work on the quark-meson model within the Functional RG suggests that the transition 


temperature becomes even smaller if one includes wave function renormalizations 


In 


spite of this, the slope of the function T c (m is roughly the same as in our study. This is an 
additional hint that neglecting the gauge degrees of freedom could indeed be responsible for 
the difference in the results from the quark-meson model compared to lattice calculations. 

A recent study in terms of the Functional RG (FRG), which incorporates gluonic degrees 
of freedom and four-fermion interactions, shows reasonable agreement with results from 
lattice studies of the chiral phase transition temperature for two and three massless quark- 
flavors (aJ. 

Results for T c in the chiral limit from various lattice and RG approaches are summarized 
in Tab. EH As can be seen from the table, there is some uncertainty in the value of the 
chiral phase transition temperature in lattice calculations, which is mainly due to different 
implementations of the fermions. 
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FIG. 2: Chiral phase transition temperature T c in infinite volume as a function of m n (T = 0). The 
dots show the result of our RG-calculation, and the line shows the result of the fit function defined 
in Eq. (I2U1) . 

IV. CHIRAL PHASE TRANSITION TEMPERATURE IN FINITE SPATIAL 

VOLUMES 

We now turn to the investigation of the chiral phase transition temperature in finite 
spatial volumes. As in section lllll we define the phase transition temperature T c via the 
minimum of the cr-mass. Putting the system in a finite volume introduces an additional 
scale. Let us first discuss the influence of this additional length scale L on the cr-mass and 
7r-mass. In Fig. |21 we show the cr-mass and 7r-mass for mi°' ) = 100 MeV and with periodic 
boundary conditions for the quarks, as a function of the temperature, for both a small 
volume L = 1 fm and a large volume L = 4 fin. The minimum of the cr-mass is clearly 
visible in the plot. Above the transition temperature T c , where chiral symmetry is restored, 
the cr- and 7r-mass are degenerate, independent of the size of the volume. In order to gain 
a better understanding of the meson masses and their dependence on the scales L and T 
in this temperature regime, we peform a perturbative one-loop calculation. Since the chiral 
phase is a non-perturbative phenomena, such a one-loop calculation is not trustworthy in the 
vicinity of the critical temperature: The extraction of the critical temperature fails, leading 
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Reference 

Method 

T c [MeV] 

this work 

Proper-time RG 

148 

Berges (1997) 

m 

Functional RG, quark-meson model 

100.7 

Schaefer (1999) 

[45] 

Proper-time RG, quark-meson model 

149 

Braun (2003) 

EH 

Proper-time RG, quark-meson model 

154 

Schaefer (2004) 

[8] 

Proper-time RG, quark-meson model 

142 

Braun (2005) 

m 

Functional RG, QCD (see caption) 

186 

Gottlieb (1996) 

E&] 

lattice 16 3 x 8 (staggered) 

128 + 9 

Karsch (2000) 

[13] 

lattice 16 3 x 4 (improved staggered) 

173 + 8 

CP-PACS (2000) 

m 

lattice 16 3 x 4 (improved Wilson) 

171 + 4 

Bornyakov (2005) [60] 

lattice 16 3 x 8 (improved Wilson) 

173 + 3 


TABLE II: Chiral phase transition temperature in the chiral limit ( m^ —> 0) from different RG 
approaches for the quark-meson model and for QCD, and from lattice simulations. We have 
restricted our choice of lattice references to the case of Nf = 2 flavors that we have treated here. 

nnn 

More recent lattice results have been obtained for Nf = 2 + 1 flavors, see e.g. .2, [lj, l&a]. The 


difference in the RG results arises from a weak dependence of T r on the initial values at the UV 


scale and on the choice of the cutoff-function Eq. (gj). In Ref. [36(, T c is calculated from a study 
which incorporates the running QCD coupling and all possible four-fermion interactions in the 
chiral symmetric regime. 


to an unphysical complex temperature |61|]. Here, higher-loop terms contribute significantly, 


as pointed out by earlier RG flow studies, eg. Ref. 


. 3 


45 


57j, and a study in terms of 


many-body resummation techniques [62|. Moreover, we neglect the quark contributions 
in this calculation, since they are suppressed by the appearance of a thermal Matsubara 
mass. In contrast, the bosonic fields have vanishing Matsubara mass and therefore their 
contributions dominate at high temperature. Thus our starting point for the calculation of 
the mass correction is the scalar 0(N = IV|)-model with the Lagrangian 


£ = ]++) 2 + j m V 2 + 4 


(23) 


where we have introduced the 0(IV)-vector </> = (</>!, ...,0 tv), and the parameters are m = 
m uv and A = A uv- 
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The mass correction 5m 2 (T, L), which is due to finite volume and finite temperature 
effects, can be decomposed into a sum of two contributions, 8m\(T,L-^oo) and 5m 2 (T, L). 
We refer to Appendix El for details of the calculation. 

First, in the regime defined by 0 the contribution 5m\(T, L —> oo) dominates. 

One can estimate 8m\(T, L— >oo) for large temperatures and volumes as 

<5mi(T, L—> 00 ) ~ (ZV + 2)— T 2 for T —■» oo. (24) 


In this case, the meson masses depend linearly on the temperature, in agreement with the 


result from Ref. 


611. 


Second, if T and L are of the same order of magnitude, the mass correction is in essence 
given by 


6m%(T,L ) 


(1V + 2)AT 

(27r)i 


EE' 

n=-oo {Zi} 


(mi ) 2 + 4vr 2 n 2 (TL ) 2 

T 2 


K 


l 2 ((mL) 2 + 47 r 2 n 2 (TL) 2 )), 
(25) 


where K n denotes the modified Bessel-functions with index n — 1. The vector l is defined 
as l = {Zi, Z 2 , ^ 3 } and the prime indicates that the term with / = 0 is excluded from the 
summation. Note that 5m, 2 (T, L ) has a complicated dependence on T and L, but we observe 
that it scales with j, rather than with T 2 . This explains the difference between the slopes 
of the meson masses in the regime defined by A > 1 > £, anc j j n the regime defined by 
0 < ^ < 1, which can be seen in Fig. E 

In contrast, in the limit TL 1 one obtains 


8m 2 2 (T,L ) fa 3 ^ N ^ J^'-^exp \Jl 2 ({mL) 2 + An 2 n 2 (TL) 2 )^j. (26) 

^ n=—00 p;} V t 

The contributions from the non-vanishing thermal Matsubara-modes to 8m 2 (T,L) drop 
exponentially, and <5m|(T, L) becomes a linear function in the temperature T, due to the 
zeroth thermal Matsubara-mode n — 0. Therefore, for TL 3 > 1, 5m\(T,L) is a sub-leading 
correction to the meson masses, compared to the contribution 5m\(T , L). This describes the 
results for periodic quark boundary conditions well. 

Similar behavior is found for anti-periodic boundary conditions of the quark fields. How¬ 
ever, there is one essential difference between periodic and anti-periodic boundary condi¬ 
tions: In the case of anti-periodic boundary conditions, the quark fields have a finite minimal 
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T [MeV] 


FIG. 3: Sigma- and Pion-mass as functions of the temperature T, with rn^ = 100 MeV and 
periodic quark boundary conditions, for L = 1 fm and L = 4 frn. The solid black and the red/gray 
lines show the pion mass for L = lfm and L = 4fm, respectively, whereas the dotted black and 
red/gray lines show the sigma mass for L = 1 fm and for L = 4 fm. 

infrared momentum 

P™ = \, (27) 

which is illustrated in Fig. 01 hi the quark-propagator, the minimal value acts as an 
additional mass term which increases for decreasing volume sizes. The quark fields decouple 
from the RG flow as soon as the IR-cutoff scale k in Eq. (HH> drops below Consequently, 
the mesons are the only dynamical degrees of freedom in the theory for k < p™ n '- The system 
becomes equivalent to an 0(4)-model, which remains in the symmetric phase for k —> 0. 

In general, the momentum scale should be much smaller than the UV cut off A w, 

2vr 

Pp = — < Xuv ■ (28) 

If Pp becomes comparable to Auv, there are no dynamical degrees of freedom left in our 
model. 
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FIG. 4: Schematic representation of the relation between the momentum summation and the 
UV- and IR-cutoff in momentum space, for anti-periodic boundary conditions (left panel) and 
for periodic boundary conditions (right panel). The UV-cutoff is denoted by A, while k denotes 
the variable IR-cutoff of the RG scheme. The arrow indicates the direction of the RG flow from 
k = A uv to k = 0. 

We now present onr main results for the volume dependence of the chiral phase transition 
temperature. Fig. 0 contains plots of the transition temperature T c as a function of the 
volume size, for different values of the pion mass at zero temperature, m^\ and for different 
choices for the quark boundary conditions. For small, realistic pion masses, mi 0) = 100 MeV, 
already for L = 4 fin the results in the upper panel of Fig. 0 show a deviation of T c from 
its infinite volume value of about 6%, independent of the choice of boundary conditions. 
For small volume sizes L < —hr, T c is strongly affected by the choice of the boundary 

m) r 

conditions for the fermions. For anti-periodic boundary conditions, we observe that T c 
decreases strongly for small volume sizes. As already discussed above, this is because of the 
additional infrared cutoff p™ in for the momenta of the quark fields. It is due to this additional 
IR-cutoff for the quarks that the system remains in the symmetric phase for small volume 
sizes (3, Q] • 

For periodic quark boundary conditions, we observe a weaker volume dependence of T c , 
since the condensation of the quarks is not prevented by the additional IR-cutoff p™ in . For 
m i° = 100 MeV and L < 1.5 fin, we observe that T c is almost independent of L, which may 
be due to the fact that p £ approaches A uv . 

For large pion masses, m ® > 300 MeV, and for L > 2 fm, T c depends only weakly on the 
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FIG. 5: The chiral phase transition temperature T c as a function of the volume size L, for different 
pion masses m^\ and for both choices of spatial quark boundary conditions. Results for periodic 
(anti-periodic) boundary conditions are denoted by circles (triangles), values for the pion mass at 
T = 0 in infinite volume are given in the panels. 
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FIG. 6: Chiral phase transition temperature T c as a function of m^(T = 0) for different box sizes. 
We show the results for L = 1.5, 2.5, 3.5 fm from bottom to top. For comparison, we also show 
the results for infinite volume (solid line) from Fig. [21 

box size. The deviation from its infinite volume value is less than 1% already for L « 2.5 fm. 
We observe only a weak dependence on the choice of the fermionic boundary conditions, as 
well. The reason is that the length scale set by the pion mass, L v ~ —for, is much smaller 
than the box size L. Therefore the volume dependence of T c is governed by the pion mass 
scale rrin\ rather than the scale set by the spatial box size: pion fluctuations are more 
strongly suppressed by their large mass than by the long-wavelength cutoff from the finite 
volume. This observation implies that lattice results for T c are not affected by the finite 
volume to any considerable degree, provided the pion mass is large, m i° > 300 MeV. 

Finally, we stress that finite volumes make the coefficient cq in Eq. © bigger for smaller 
volumes. This can be seen from Tab. cm and Fig. [HI where the slope of T c (rri£\ L) as a 
function of is even larger at smaller values of L. 


V. CONCLUSIONS 


In this article, we have presented results for the volume dependence of the phase transition 
temperature for the chiral symmetry restoration transition. Our investigation is based on 
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L [fm] 

m = 100 MeV 

m 1°^ = 200 MeV 

m\ r°^ = 300 MeV 

1.5 

79.0 MeV 

164.3 MeV 

225.7 MeV 

2.5 

111.3 MeV 

197.1 MeV 

246.6 MeV 

3.5 

157.9 MeV 

206.3 MeV 

249 MeV 

00 

178.1 MeV 

208.3 MeV 

249.3 MeV 


TABLE III: Chiral phase transition temperatures T c (mi°\ L) as function of the pion mass m ® 
and the box size L for periodic boundary conditions for the quark-fields. For comparison, the 
corresponding values in infinite spatial volume are also given. 


the quark-meson model and uses Renormalization Group methods. In this way, we obtain 
non-perturbative results for the transition temperature for various values of the pion mass. 

In general, no phase transition can occur in a system of finite volume. The evaluation of 
lattice results therefore uses a scaling analysis in quark mass and temperature jl^ . [ 47 I . Idil Ifiiil] . 


and recently finite-size scaling jl^j] as an analytical tool. We expect that a Renormalization 
Group analysis of the critical behavior can complement these approaches. 

In the present paper, we have focused on the chiral phase transition temperature, which is 
not universal and also model-dependent. However, the relative shift of the temperature from 
infinite to finite volume should depend mainly on the pion mass and the pion decay constant, 
which are independent of the model and represent an external input to our calculations. 

We find that finite volume effects for the transition temperature remain small for large 
pion masses m n > 300 MeV, as long as the volume is of the order L > 2 fm in the spatial 
directions. The scale for the appearance of sizable finite-volume effects is given by the pion 
mass r , and the effects remains small as long as L 

On the other hand, finite size effects are sizable already at a lattice extent of L ~ 2 fm for 
realistic pion masses of the order of 100 MeV. We expect therefore that finite volume effects 
will become more relevant in future simulations with realistic pion masses. The strategy 
for lattice calculations should then be to simulate in volumes where the value m n L is large 
enough to keep the finite volume effects down to an acceptable size, and to extrapolate to 
smaller pion masses and the chiral limit. 

We note that the choice of periodic boundary conditions in spatial directions for the 
quark fields, which is commonly employed in lattice simulations, leads to a much smaller 
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finite volume effect on the transition temperature. This conclusion a gree s with our results 
for the volume dependence of the pion mass and pion decay constant 26|. 

The dependence of the transition temperature on the quark mass in the quark-meson 
model is much stronger than that observed in lattice simulations. This must remain a 
puzzle in the present approach and calls for a more extensive consideration of the gluon 


136 ]. 


dynamics. Work in this direction has already been started 1361 . In concert with other non- 
perturbative methods, the Renormalization Group approach should prove to be valuable as 
a complement to QCD lattice calculations. 
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APPENDIX A: INITIAL VALUES AT THE UV-SCALE 


In this section, we give an overview of the initial values for the flow equations of the 
coefficients dij(k). As already discussed in Sec. |TT1 we fix the values of the coefficients 
dij(Auv) at the scale A jjy in the infinite four-dimensional Euclidean volume in such a way 
that the values of m^ and f n given by chiral perturbation 123] are reproduced. However, some 
freedom remains in the choice of the starting values for the coefficients. But as discussed in 

n 

26|, different sets of parameters give the same dependence on the size of the four dimensional 
volume, provided that they lead to the same values of the pion decay constant and pion mass 
in infinite volume. The values which we have used for our calculation can be determined 
from Tab. El In this table, we give the value for the parameter muy. Together with 
the parameter A uv — A(A[/y) and the evolving minimum of the potential <7o(/c), the two 
coefficients that are initially non-zero can be expressed as 


ooi(Af/v) — -{m 2 uv -\-\uv^l{A uv )) 

O02(Af jy) = -A uv ■ 


(Al) 

(A2) 
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muv [MeV] 

gm c [MeV] 

m 7r (k —> 0) [MeV] 

f n (k - 0) [MeV] 

794 

0.036 

30 

87 

792 

0.39 

50 

87 

788 

1.08 

75 

88 

779 

2.10 

100 

90 

775 

3.43 

125 

91 

767 

5.15 

150 

93 

757 

7.275 

175 

95 

748 

9.85 

200 

97 

737 

12.93 

225 

99 

725 

17.00 

250 

101 

711 

20.80 

275 

103 

698 

25.70 

300 

105 


TABLE IV: Boundary conditions aij(Ajjv) for the differential equations for the coefficients a,ij(k). 
The corresponding values for m, and which result from the RG flow at zero temperature and in 
infinite volume are shown in the last two columns. The coefficients aoi and ao 2 ) which are initially 
non-zero, can be expressed in terms of the two parameters mjjv and A;yy, for the relation see text. 
As an initial condition for the flow equation for A (k), we have chosen A(A uv) = A uv = 60 for all 
parameter sets. The initial conditions for the differential equations for the remaining coefficients 
a,ij{k ) of the potential Eq. (ffifll are ciij(Auv) = 0. In addition, we use g = 3.258 for the Yukawa- 
coupling. In our notation, gm c is the current quark mass. 

All other coefficients are initially set to zero. In principle, one should choose <7o(A-uv) = 0 
as initial condition, but in order to avoid numerical problems at the UV-scale, we choose 
a small but finite value for <r 0 at the UV-scale, e.g. a 0 (A uv ) = 0.1 MeV. The parameters 
muv and A uv can be interpreted as the usual meson masses and four-point, couplings at the 
UV scale. The meson masses are related to the coefficients a,-,- and the minimum <7 n of the 

n 

potential by |25|| 

ml = 2a 0 i 

ml = ml + 2a 2 o + 4an<To + 8ao 2 o'o 


(A3) 

(A4) 
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APPENDIX B: ONE-LOOP CALCULATION OF THE MESON MASSES 

In order to gain a better understanding of the slope of the meson masses as a function of 
temperature in the symmetric phase, we compared the RG results to a one-loop calculation 
for the mass of a scalar 0(A’)-model, as defined by the Lagrangian in Eq. f!2.‘lll . More details 
of the calculation are given here. 

The effective potential in one-loop approximation in a d-dimensional Euclidean space 
reads 

rlL/ 


U 1L (0) = -Trlog (—«9 2 + M 2 (0)) 


N 


1 T v-x- ds 

'2^^^ ■ 7 exp 

3= 1 {n<} J A? 


+ + (Bl) 

i= 1 


where we have used Schwinger’s proper time representation of the logarithm, supplemented 
by a UV-cutoff A. In addition, we have introduced the abbreviations M 2 (</>) = m 2 + 3A</> 2 and 
Mj (0) = m 2 + A0 2 for j = 2,..., N. The sum runs over all permutations of {ni,..., n^}. Here 
L denotes the length of the d -dimensional box in space directions, and A denotes the length 
of the box in the Euclidean time direction. In order to calculate the mass of the scalar field, 
we use the Jacobi-Elliptic-Theta function $ p (s), defined in Eq. (HI. Both representations of 
this function given in Eq. m deserve some comments: The first representation in Eq. HU) 
is essentially the usual Matsubara summation, where a truncation of the Matsubara sum can 
be used to perform a high-temperature or small-volume expansion of the one-loop effective 
potential. In order to get the second representation in Eq. HUD, we have applied Poisson’s 
formula to the first representation. This second representation allows to separate the UV 
divergence of the effective potential: the divergent infinite-volume contribution, which is 
given by the first term of the left-hand side of Eq. nsD, can be separated from the volume 
and temperature dependent parts. In addition, a truncation of the sum in the second 
representation can be used to perform a low-temperature or large-volume expansion of the 
one-loop effective potential. 

With the second representation of the Jacobi-Elliptic-Theta function from Eq. HSD , we 
can divide the one-loop effective potential in three contributions, 


U 1L = U£(T^O,L^oc) + UI l (T,L^oc) + U} l (T,L) 


t 1L 


1L / 


T 1L / 


(B2) 


We do not consider the contribution in the following, since we are only interested in the 
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finite-temperature and finite-volume corrections to the mass. We start with the calculation 
of the mass correction due to the volume-independent contribution Z/, 1 L (T) = f/[ 1 L (T, L —> cxj) 
which is given by 



We stress that the regulator A of the Schwinger proper-time integral can be removed here, 
A —> oo. In order to calculate the mass correction in the symmetric phase, we have to 
take the second derivative of U\ h (T) with respect to the field 0 and evaluate the resulting 
expression at </> = 0 : 


d 2 


6mi{T) = w ui {T) 


2 (N + 2)A ^<*=2 

- (m 


b =0 


(2vr; 


=i <? 2 “ J 


(B4) 


Now we make use of the integral representation of the modified Bessel Functions, 


/T\ U r 00 - 7 T 2 

2(-) AA(x) = j f . 

For T —> cx) and d = 4, we use 


(B5) 


A’,(x) 



(i « 1), 


(B6) 


and obtain the simple expression 


8m\(T) 


(N + 2)A t2 
12 


(B7) 


fill. 


for the correction 5m, 2 (T), which agrees with the result found in 

Now we turn to the calculation of the mass correction due to the contribution [/ 2 1 L (T, L). 
Since we are interested in small spatial volumes and high temperatures, we use the Poisson- 
representation for the spatial contributions, but we still use the usual Matsubara sum for 
the thermal contribution. The contribution f/ 2 L (T, L) to the potential is 


I/ 2 1 l (T,L) = 


T 


N 


2(4tt 


d— 
\ 2 




j =1 n =-oo {l t } 


ds 

P 

S 2 


s(M/(0)+4vrVT i 


f 2 L 2 

4s 


(B 8 ) 


The vector l is defined as l = {Zi, Z 2 , Z 3 }, and the prime indicates that the term with Z = 0 is 
excluded from the summation. The regulator A can also be removed here. Taking the second 
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derivative with respect to the fields <fi and evaluating at 0 = 0 , we obtain the corresponding 
mass correction 


5m 1 2 3 4 5 6 7 8 9 10 11 2 (T,L ) = 


*-U 

d^ 2 


1L 


(T,L) 


(B9) 


t >=0 


(iV + 2)ATL-V 


(2vr 


d -1 
\ 2 


EE' 

n=-oo 


m 2 + 47r 2 n 2 T 2 

l 2 


d -3 
4 


K d -3 
2 


l 2 ((mL) 2 + 47 r 2 n 2 (TL) 2 ) 


We used again the integral representation of the modified Bessel functions Eq. (HU). Finally, 
we show that this contribution becomes proportional to the temperature j for TL > 1. For 
this purpose, we use the asymptotic expansion of the Bessel-functions for large arguments, 
given by 


KJx) « W-e x (x > 1). (BIO) 

V x 

Using Eq. Em the mass correction Sm 2 (T , L) for EL 1 and d — 4 reads 


5m 2 2 (T,L) 


3(N + 2)A T 
8^7r L 


E E'yf ex p 

i=-oo pi} V l 


l 2 ({mL) 2 + 4n 2 n 2 (TL) 2 


(BH) 


which drops exponentially in the limit TL —>• cx) for non-vanishing thermal Matsubara- 
modes, whereas the contribution from the zeroth thermal Matsubara-mode (n = 0) remains 
finite and is proportional to the temperature T. For Th > 1, the mass-correction 5m 2 (T, L) 
is therefore sub-leading, compared to the contribution 5m\(T). 
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